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Abstract 

Consider the random sequential packing model with infinite input and in 
any dimension. When the input consists of non-zero volume convex solids 
we show that the total number of solids accepted over cubes of volume A is 
asymptotically normal as A — > oo. We provide a rate of approximation to 
the normal and show that the finite dimensional distributions of the packing 
measures converge to those of a mean zero generalized Gaussian field. The 
method of proof involves showing that the collection of accepted solids satisfies 
the weak spatial dependence condition known as stabilization. 



1 Main results 

Given d G N and A > 1, let L r 1A) U 2) \, ... be a sequence of independent random d- 
vectors uniformly distributed on the cube Q\ := [0, \ l l d ) d . Let S be a fixed bounded 
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closed convex set in M. d with non-empty interior (i.e., a 'solid') with centroid at the 
origin of M. d (for example, the unit ball), and for ieN, let S i: \ be the translate of 
S with centroid at U^\. So S\ := (Si t \)i>i is an infinite sequence of solids arriving 
at uniform random positions in Q\ (the centroids lie in Q\ but the solids themselves 
need not lie wholly inside Q\). 

Let the first solid S lt \ be packed, and recursively for % — 2, 3, . . ., let the i-th solid 
Si, a be packed if it does not overlap any solid in {S^\, . . . , Si-i t \} which has already 
been packed. If not packed, the i-th solid is discarded; we sometimes use accepted 
as a synonym for 'packed'. This process, known as random sequential adsorption 
(RSA) with infinite input, is irreversible and terminates when it is not possible to 
accept additional solids. At termination, we say that the sequence of solids S\ jams 
Q\ or saturates Q\. The jamming number N\ := N\(S\) denotes the number of 
solids accepted in Q\ at termination. We use the words 'jamming' and 'saturation' 
interchangeably in this paper. 

Jamming numbers N\ arise naturally in the physical, chemical, and biological 
sciences. They are considered in the description of the irreversible deposition of 
colloidal particles on a substrate (see the survey [1] and the special volume [20]), 
hard core interactions (see the survey [7]; also [25]), adsorption modelling (see [3] 
and the survey [24]) and also in the modelling of communication and reservation 
protocols (see [4, 5]). 

The extensive body of experimental results related to the large scale behavior 
of packing numbers stands in sharp contrast with the limited collection of rigorous 
mathematical results, especially in d > 2. The main obstacle to a rigorous math- 
ematical treatment of the packing process is that the short range interactions of 
arriving particles create long range spatial dependence, thus turning N\ into a sum 
of spatially correlated random variables. 

In the case where d — 1 and S — [0,1], a famous result of Renyi [21] shows that 
jamming limit, defined as lim^oo A _1 EiV A , exists as an integral which evaluates to 
roughly 0.748; also in this case, Mackenzie [10] shows that liniA^oo A _1 VariV.x exists 
as an integral which evaluates to roughly 0.03815. Dvoretzky and Robbins [6] show 
that the jamming numbers N\ are asymptotically normal as A — > oo, but their 
techniques do not address the case d > 1. 

Since the above results were established in the 1960s, progress in extending them 
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rigorously to higher dimensions has been slow until recently. Penrose [11] establishes 
the existence of a jamming limit for any d > 1 and any choice of S, and also [12] 
obtains a CLT for a related model (monolayer ballistic deposition with a rolling 
mechanism) but comments in [12] that 'Except in the case d — 1 ... a CLT for 
infinite-input continuum RSA remains elusive.' 

In the present work we show for any d and S that A^VarA^ converges to a 
positive limit and that N\ satisfies a central limit theorem, i.e., the fluctuations of 
the random variable N\ are indeed Gaussian in the large A limit. This puts the recent 
experimental results and Monte Carlo simulations of Quintanilla and Torquato [22] 
and Torquato (ch. 11.4 of [25]) on rigorous footing. We also provide a bound on the 
rate of convergence to the normal, and on the rate of convergence of A _1 E A^ A to the 
jamming limit. 

Throughout Af(Q, 1) denotes a mean zero normal random variable with variance 
one. 

Theorem 1.1 Let S\ be as above and put N\ := N\(S\). There are constants 
fj, := fi(S, d) G (0, oo) and a 2 := a 2 (S, d) G (0, oo) such that as A — > oo we have 

[X^ENx- fi\ = 0{X' 1/d ) (1.1) 

and A _1 VarA^A — > a 2 with 

= 0((\og\) 3d \~ 1/2 ). (1.2) 



sup 



N\ — E N\ 
(VarAT A )V2 



< * 



- P[Af(0, 1) < f] 



The process of accepted solids in Q\ induces a natural random point measure v\ 
on [0, l] d given by 

oo 

v\ := $^^ A_1/d ^,A 1 {5 l , A is accepted} ( L3 ) 

i=l 

where S x stands for the unit point mass at x. It also induces a natural random 
volume measure v' x on M d , normalized to have the same total measure as v\, defined 
for all Borel A C M. d by 

v' x (A) := A | A n (|J[^ _1/ "^,a :i > 1, S hX is accepted]) | (1.4) 
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where | • | denotes Lebesgue measure and A x l d A := {A 1 l d . 



x : x E A}. The measure 



u' x is not necessarily supported by Qi due to boundary effects, but for A > 1 it is 
supported by Q^, where we set := [— l,2) d (a fattened version of Qi). 

Let V\ := v\— E[z/ A ] and v' x := v' x — E[i/jJ. Let lZ(Qf) denote the class of bounded, 
almost everywhere continuous functions on Q^. For / E lZ(Qf) and /i a signed 
measure on R d with finite total mass, let (/,//) := J R fdfx. The following theorem 
provides the limit theory (law of large numbers and central limit theorems) for the 
integrals of test functions / E 1Z(Qf) against the random point measure v\ and the 
random volume measure v' x induced by the packing process. In particular, it shows 
that the finite dimensional distributions of the centered packing point measures (^a)a 
converge to those of a certain mean zero generalized Gaussian field, namely white 
noise on Qi with variance a 2 per unit volume, and likewise for the centered packing 
volume measures {p' x )\. 

Theorem 1.2 Let ji and a 2 be as in Theorem 1.1. Then for any f,g in 1Z{Qi ), 



Also, the finite- dimensional distributions of the random field (A _1//2 (/, ux), f E 1Z(Qi)) 
converge as A — > oo to those of a mean zero generalized Gaussian field with covari- 
ance kernel 



Moreover, the same conclusions hold with v\ and v\ replaced by u' x and v' x respec- 
tively. 

Remarks. 

1. Finite input. Let r E (0, oo) and let \x] denote the smallest integer greater 
than or equal to x. Inputting only the first [At] solids of the sequence S\ yields 
RSA packing of the cube Q\ with finite input. The finite-input packing number, 
i.e., the total number of solids accepted from S\ t \, S^a, S^tA^a, is asymptotically 




and 
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normal as A — > oo with r fixed. This is proved in [17], and extended in [2] to the 
case where the spatial coordinates come from a non-homogeneous point process. 
Packing measures induced by RSA packing with finite input have finite dimensional 
distributions converging to those of a mean zero generalized Gaussian field with a 
covariance structure depending upon the underlying density of points [2]. 

2. Stabilization. One might expect that the restriction of the packing measure i>\ 
or v' x to a localized region of space depends only on incoming particles with 'nearby' 
spatial locations, in some well-defined sense. This local dependency property is 
denoted stabilization; when the region of spatial dependency has a diameter with 
an exponentially decaying tail, it is called exponential stabilization. These notions 
are spelt out in general terms in Section 2. Theorem 2.1 provides a general spatial 
limit theory for exponentially stabilizing measures; this is an infinite-input analog 
to known results [2, 13, 14, 15] for the finite-input setting, and is of independent 
interest. 

A form of stabilization for infinite input RSA was proved in [11], but without 
any tail bounds. Exponential stabilization in the infinite input setting is perhaps not 
surprising, but it has been challenging to rigorously establish this key localization 
feature. In Section 3, we show that infinite-input packing measures stabilize expo- 
nentially, so that the general results of Section 2 are applicable to these measures. 

3. Related models in the literature (see e.g. [17]) include cooperative sequential 
adsorption, RSA with solids of random size or shape, ballistic deposition with a 
rolling mechanism, and spatial birth-growth models. For all of these models, limit 
theorems in the finite-input setting are discussed in [12]. It seems likely that these 
can be extended to the infinite-input setting using the methods of this paper, al- 
though we do not discuss any of them in detail. Nor do we consider non-homogeneous 
point processes as input. 

4. Rates of convergence. Even in d — 1, the rate given by Theorem 1.1 is new. 
Quintanilla and Torquato [22] use Monte Carlo simulations to predict convergence 
of the distribution function for N\ to that of a normal, but they do not obtain 
rates. Penrose and Yukich [19] obtain rates of approximation to the normal for RSA 
packing with finite (Poisson) input. 

5. Numerical values. We do not provide any new analytical methods for com- 
puting numerical values of /i and a 2 when d > 1. 
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6. Jamming variability. A significant amount of work is needed (see Section 4) 
to show that the limiting variance a 2 in Theorems 1.1 and 1.2 is non-zero, and we 
prove this using the following notions. 

Given L > 0, we shall us say that a point set rj C R d \ [0, L] d is admissible if the 
translates of S centered at the points of rj are non-overlapping. Given such an rj, let 
N[[0, L] d \r/] denote the (random) number of solids from the sequence S L d which are 
packed in [0, L] d given the pre-packed configuration rj. In other words, N[[0, L] d \rf\ 
arises as the number of solids packed in [0, L] d in the course of the usual infinite input 
packing process subject to the additional rule that an incoming solid is discarded 
should it overlap any solid centered at a point of rj. Say that the convex body S has 
jamming variability if there exists a L > such that inf,, VariV[[0, L] d \rf\ > with 
the infimum taken over admissible point sets i] C M. d \ [0, L] d . 

In Proposition 4.1 we shall show that each bounded convex body S C M. d with 
non-empty interior has jamming variability. 

7. We let ds stand for the diameter of S. In our proofs, we shall assume that 
2ds < 1. This assumption entails no loss of generality, since once we have proved 
Theorems 1.1 and 1.2 under this assumption, the results follow for general S by 
obvious scaling arguments. 

2 Terminology, auxiliary results 

Let M + := [0, oo). Given a point (xi, . . . , x d , t) = (x,t) G R d x R + , the first d 
coordinates of the point will be interpreted as spatial components with the (d+ 1)- 
st regarded as a time mark. Let us say a point set X C M. d x M + is temporally locally 
finite (or TLF for short) if XC\ (M. d x [0, t]) is finite for all t > 0. Loosely speaking, X 
is TLF if it is finite in the spatial directions and locally finite in the time direction. 

In this section we adapt the general results and terminology from [2, 14, 15, 19] 
on limit theory for stabilizing spatial measures defined in terms of finite point sets 
in R d , to to the setting of spatial measures defined in terms of TLF point sets 
in M, d x R + (typically obtained as Poisson processes). In subsequent sections, we 
show that these general results can be applied to obtain the limit theorems for RSA 
described in Section I. 

For x E M. d and r > 0, let B r (x) denote the Euclidean ball centered at x of 
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radius r. We abbreviate B r (0) by B r . Given X C R d x R+, a > and y G M d , we 
let y + a^Y := {(y + ax, t) : (x,t) G X}; in other words, scalar multiplication and 
translation on M. d x M + act only on the spatial components. For A C M. d we write 
y + a A for + ax : x G A}; also, we write cM for the boundary of A, and write A + 
for A x R + . For nonempty subsets A, A' of M d , write D 2 (A,A') for the Euclidean 
distance between them, i.e. D 2 (A, A') := inf{|x — y\ : x <E A,y e A'}. 

Let be an M^-valued function defined for all pairs (X, A), where X is a 

TLF subset of W 1 x M + and A is a Borel subset of W 1 . Throughout this section we 
make the following assumptions on £: 

1. £(-,A) is measurable for each Borel A, 

2. £(X, •) is a finite measure on M d for each TLF ^cK d xR + , 

3. £ is translation invariant, that is £(« + Af, i + A) = A) for all i G Z d , all 
TLF X C M d x R + , and all Borel A C R d , 

4. ^ is uniformly locally bounded (or just boundediox short) in the sense that there 
is a finite constant ||£||oo sucn tnat f° r & H TLF A 1 C l d x R + we have 

i{X,%\f)<\\i\U. (2.1) 

5. £ is locally supported, i.e. there exists a constant p such that £(Af, A) = 
whenever D 2 (X,A) > p. 

Note that if £(X, ■) is a point measure supported by the points of X, then £ is locally 
supported (in fact, in this case we can set p = 0). 

For all A > 0, let V\ denote a homogeneous Poisson point process in M d x M + 
with intensity measure \dx x ds, with <ir denoting Lebesgue measure on M. d and cLs 
Lebesgue measure on R + . We put V :—V\. 

Thermodynamic limits and central limit theorems for functionals in geometric 
probability are often proved by showing that the functionals satisfy a type of local 
spatial dependence known as stabilization [2, 13, 14, 15, 17, 18, 23] and that will be 
our goal here as well. First, we adapt the definitions in [2, 13, 14] to the context of 
measures defined in terms of TLF point sets in M. d . Recall that Q\ denotes the cube 
[0,A 1 / d ) d . 
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Definition 2.1 We say £ is homogeneously stabilizing if there exists an a.s. finite 
random variable R' (a radius of homogeneous stabilization for such that for all 
TLF X C (lR d \ B R ,)+ we have 

a(v n (b r ,) + ) u x , go = £{v n (b r ,) + , q 1 ). (2.2) 

VFe say £ is exponentially stabilizing if (i) it is homogeneously stabilizing and R' 
can be chosen so that limsup i _ +00 L' 1 log P[R' > L] < 0, and (ii) for all A > 1 and 
all i G Z d ; i/iere eaxste a random variable R := R^(i, \) (a radius of stabilization for 
£ at i with respect to V in (Q\) + ) such that for all TLF X C [Q\ \ B R {i)] + , and 
all Borel AC. Qi, we have 

a(vn [B R {i) n Q x ] + ) u x,i + A) = £ (v n n Q A ]+, i + A) (2.3) 

and moreover the tail probability r(L) defined for L > by 

r(L):= sup P[R*(i,\)>L] (2.4) 

A>1, ieZ d 

satisfies limsup^^ logr(L) < 0. 

Loosely speaking, R := R^(i,X) is a radius of stabilization if the ^-measure on 
i + Qi is unaffected by changes to the Poisson points outside B R (i) (but inside Q\). 
When £ is homogeneously stabilizing, the limit 

£(7>,i + Qi) := lim£(7>n (£,.(;))+, i + Qi) 

exists almost surely for all % G 2, d . The random variables (£(P, % + Qi),i G Z d ) form 
a stationary random field. 

Given £, for all A > 0, all TLF C R d x R+, and all Borel AcK d we let 
£x(X, A) := £(\ 1/d X, \ 1/d A). Define the random measure fi{ on R d by 

£(-)--=WP*riQ ir ) (2.5) 

and the centered version jfi x := ^ — E[/4]. By the assumed locally supported 
property of £, /i A is supported by the fattened cube := [—1, 2) d for large enough 
A. 
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If £ is stabilizing, define //(£) := E [£(P, Qi)] and and if £ is exponentially stabi- 
lizing, define 

^(O^^Cov^QO, £(P,i + Qi)], 

where the sum can be shown to converge absolutely by exponential stabilization 
and (2.1). The following general theorem provides laws of large numbers and normal 
approximation results for (/, fi\), suitably scaled and centered, for / G lZ(Qf). This 
set of results for measures determined by TLF point sets is similar to previously 
known results for measures determined by finite point sets (Theorem 2.1 of [18], 
Theorem 2.1 of [2], Theorem 2.3 of [2], and Corollary 2.4 of [19]). 

Theorem 2.1 Suppose that £ is exponentially stabilizing. Then as A — > oo ; for f 
and g in TZ.(Qf) we have 



and 



Also, 



lim X-'E [(/, /4>] = //(£) / f(x)dx 
hm X^CavKf,^, (</,/4>] = o- 2 (0 [ f(x)g(x)dx. 



|A- 1 E[ / xi(g+)]-MOI = o(A- 1 / d ) 

Moreover, ifcr 2 (t;) > then 

"/4(Q+)-E[/4(Q+)] 



(2.6) 

(2.7) 
(2.8) 



sup 



P 



(Var[/xi(g+)])V2 



< * 



iW(o, i) < f] 



= 0((logA) M A- 1/2 ) (2.9) 



and the finite- dimensional distributions of the random field (A -1 / 2 (/, /^), / G 7£(<3i~)) 
converge as A — > oo to inose o/ a mean zero generalized Gaussian field with covari- 
ance kernel 

(/, 9) ^ ^(0 / f(x)g(x)dx, /, 5 G rc(Q+). 



We shall use Theorem 2.1 to prove the results on RSA described in Section 1. 
It seems likely that Theorem 2.1 can also be applied to obtain similar results for 
the related models listed in Remark 3 of Section 1. For some of these, certain 
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generalizations of Theorem 2.1 may be needed; for example, in some cases one may 
need to allow for the Poisson points to carry independent identically distributed 
random marks, and in others the boundedness condition (2.1) may need to be relaxed 
to a moments condition. It seems likely that little change to the proof of Theorem 
2.1 will be needed to cover these generalizations. 

As we shall see shortly, the thermodynamic limits (2.6) and (2.8) do not require 
exponential decay of the stabilization radius for £, but in fact hold under weaker 
decay conditions. We expect that (2.7) also holds under weaker decay conditions 
on the stabilization radius, and also that the boundedness condition (2.1) can be 
relaxed to a moments condition in Theorem 2.1, but for simplicity we shall assume 
throughout that £ is exponentially stabilizing and satisfies (2.1). Also, if we restrict 
attention to / supported by Qi, we do not need the condition that £ be locally 
supported. 

The rest of this section is devoted to proving Theorem 2.1. We shall use the 
following notation. Given / G 1Z(Qf), we extend / to the whole of R d by setting 
f(x) =0foriel d \ Q+. Given TLF X C (M. d ) + , and A > 0, write (/,&(*)) 
for f Rd f (x)£\(X , dx) (the integral of / with respect to the measure £ X (X,-)). For 
j G \- l,d Z d , let f XJ : R d -> R be given by f XJ (x) = f(x) for x E j + Q 1/x , and 
fxj(x) = otherwise. Then 



f(X,j) := sup{f(x) :xEj + Q 1/x }; f(X,j) := m£{f(x) :xEj + Q 1/x }. 
For x e R d let i\(x) be the choice of % G \~ 1/d Z d such that x G i + Qy x . 

Proof of (2.6). Let / G K(Qt). Then by (2.10), we have 




(2.10) 



je\- 1 / d i d 



Also, let 



A- 1 E[(/,^>] = A- 1 n(fx„UVxn(Q 1 ) + ))] 




(2.11) 
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For x G W 1 \ dQi, with / continuous at x, we assert that as A — > oo, 

E KfxMzhWPx n (Q0 + )>] - |i(0/(x)l gi (x). (2.12) 

This clearly holds for x G lR d \ [0, l] d , since both sides are zero for large A, by the 
locally supported property of £. To see (2.12) for x G (0, l) d , observe that the left 
side has the upper bound 

E [{fxMxhWPx H {Qi)+))] < 7(A, i A (x))E [6 n (Q0+, i A (x) + Q 1/A )] 

= J(A, ix(x))E [£(V n (Q A )+, hiX^x) + (2.13) 

and has a similar lower bound with f(X,i\(x)) instead of f(\,i x (x)). If / is con- 
tinuous at x, then both f(X,i\(x)) and f(\,i\(x)) tend to f(x), so to prove (2.12) 
it suffices to show the expectation in the last line of (2.13) converges to //(£)• By 
translation invariance, this expectation equals 

E[aVn(-i 1 (X 1 / d x) + Q,) +: Q 1 )}. 

For x in the interior of Q±, the set —ii(X 1 ^ d x) + Q\ has limit set K d as A — > oo, i.e. 
for any r < oo the ball £? r is contained in — ii{X}l d x) + Q\ for large enough A. Hence 
by stabilization, 

Z(v n HiCA 1 '"*) + Q A )+, Qi) QO (2.14) 

and by (2.1), the corresponding expectations converge. This demonstrates (2.12). 

The integrand in (2.11) is dominated by a constant for x G Qf, and is zero for 
x ^ Qf. So by (2.12) and dominated convergence applied to (2.11), we obtain (2.6). 
□ 

Proof of (2.8). For this proof, set f(x) = 1 on . We need to bound the error 
term in (2.6) for this choice of /, which we do by using (2.11) again. For x G M d , 
let X(x, A) be the integrand in (2.11), i.e. set X(x, A) := (f\ Mx) , £\{V\ n (Qi)+)) 
with our current choice of /; also set Y{x, A) := £(V fl (— ii{X 1 / d x) + Q\) + , Qi). If 
x G (0, 1 - X-^y, then 

E [X(x, A)] = E [a (Pa n (Qi)+, i A (x) + Q 1/x )\ 
= E[e(?n(Q A ) +)Jl (A 1/d x) + gi)] 
= E UV n (-^(A 1 /^) + Q A ) + , Q 1 )} = E [y(x, A)]. (2.15) 
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Abbreviating the Euclidean distance D 2 ({y},A) by D 2 (y,A) is we have 

D 2 (0,d(-i 1 (X 1 / d x) + Q x )) = D 2 {ii(\ l/d x), dQ\) 
> D 2 (X 1/d x, dQx) -Vd = X 1/d D 2 (x, dQ ± ) - Vd. 

Hence the ball B xl/dD ^ x d Q^_^ is contained in the box — ^(A^rr) + Q\, so with R' 
denoting the radius of homogeneous stabilization of £, 

Y(x, \)1{R' < X 1/d D 2 (x, dQj -Vd}= £(V, Qi)l{# < X 1/d D 2 (x, dQ,) - Vd}. 

(2.16) 

Set n := = E Qi)]. By (2.15), (2.16) and (2.1), we have for x £ (0, 1 — 
\- 1 / d ) d that 

|E [X(x, A)] -fi\ = \E [Y(x, A)] - n\ 
= |E[(y(x,A) -£(V,Qi))l{B! > X 1/d D 2 (x, dQi) - Vd}}\ 
< 2U\\ 00 P[R' > X 1/d D 2 (x, dQi) - Vd\ 

and so by exponential stabilization, there is a constant K > such that 

|E [X(x, A)] - n\ < Kexp(-X 1/d D 2 (x,dQ 1 )/K). (2.17) 

Also by (2.1), for suitable K the same bound (2.17) for holds trivially for 
x e Qi \ (0,1 - X^/y, and hence (2.17) holds for all x E Qi. By (2.17), it is 
straightforward to deduce that 



/ \E[X(x,X)}- fi\dx = 0(X' 1/d ). (2.18) 
JQi 



Also, for x E R d \ Qi with D 2 (x, dQx) > X~^ d we have E [X(x, A)] = 0, and X(x, A) 
is uniformly bounded by (2.1), so that 

/ \E[X(x,X)}\dx = 0(X- 1/d ). 

JR d \Qi 

Combining this with (2.18) and using (2.11) gives us (2.8). □ 

Proof of (2.7). Let / G lZ(Qf) and assume / is nonnegative. By linearity, it 
suffices to prove (2.7) in the case where / is nonnegative and / = g, so we now 
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assume this. First, we assert that there is a constant K, independent of A, such that 
for all A > 1 and all % G A" 1 /^, z G Z d , we have 

|Cov[(/ Ai ,,6(Pn(g 1 ) + )),(/ w - 1/d2 ,6(^n(g 1 ) + ))]| <Kexp(-\ z \/K). (2.19) 

This can be proved by arguments similar to those in, e.g., the proof of Lemma 4.1 
in [2] or that of Lemma 4.2 in [15]. By (2.10), we have 

A- 1 Var[(/,^)] = A" 1 £ Cov[(/ A|i , U^x n (Qi)+)>, (/a,„6(^a n (Qi)+)>] 

i,je\- 1 / d 1 d 

dx^Cov[(/ AilA(:E) ,e(PAn(g 1 ) + )) ) (/A, A ( :c ) + A-^,e(^An(g 1 ) + )))] 

(2.20) 

where the inner sum converges absolutely by (2.19) and is zero for x Qf. 

Fix x G (0, l) d and z G Z d , with / continuous at x. Then we have the upper 
bound 

E[(f XMx) ,WPx n (Qi)+)>(/A )iA (x)+A-v^,a(^ n (Qi)+)>] 

<7(n(^),A)7(^A(^) + A- 1 /^,A) 
xE [Ca(Pa n (Q0+, i x (x) + Qi/a)6(^a n (Q0 + , i x (x) + \- l ' d z + Q 1/A )] (2.21) 

and a similar lower bound with f(i x (x),X) f(i\(x) + \~ x l d z, A) replaced by 
/(u(x), A) + A- 1 /^, A). Note that both J(i x (x), A) f(i x (x) + X^^z, A) and 

f(i\(x), A) f(i\(x) + \~ 1 / d z, A) converge as A — > oo to f 2 (x). 
By scaling and translation invariance of £, we have 

E n (Qi)+, i A (x) + Qi/a)&(7>a n (Qi)+, i A (x) + A- 1/d * + Q 1/x )\ 
= e [£(p n (Q A )+, n(A 1/d x) + Q0e(P n (Q x ) + , hiX^x) + z + Q0] 

= e [f (7? n (-^(A 1 ^) + q a ) + , goe(P n Hi(A 1A V) + q x ) + , z + Qi)]. 

By a similar argument to (2.14), as A — > oo we have 

e(pn(-i 1 (A 1 /^)+g A ) + , g 1 )e(Pn(-^ 1 (A 1 /^)+g A ) + , z+qo ^> £(p, *+Q0 

and since £ is bounded (2.1), the expectations converge. Hence, by (2.21) and the 
similar lower bound, 

E KfxMxhWP* n (Qi) + )>(/A >iA (x)+„^(^A n (Qi)+))] 

- / 2 (x)E [£(7>, gO£(P, z + Q0]. (2.22) 
13 



Also, E{(f x , h ( x ),£\(V\ n (Qi)+))] converges to f(x)fj,(Z) by (2.12), and a similar 
argument yields 

E[(/a, a (,) + A-V^,6(^A n (Qi)+)>] - 

Combining these with (2.22), we obtain that for a; G (0, l) d with / continuous at x, 

hm Cov [(/ A)iA(a! ),6(PA n (Qi)+)>, (/A, iA (, )+ „6(^A n (Q0+))] 

A— »oo 

= f 2 (x)Cov [£(V, QiU{V, z + Qi)] l 0l (x). (2.23) 

Also, (2.23) holds for x G IR d \ [0, l] d as well, since both sides are zero for large A. 
By (2.19), (2.23) and the dominated convergence theorem, applied to the last line 
of (2.20), we obtain 

hm A- 1 Var[(/,^>]=(7 2 (0 

A— »oo 

In other words, we have demonstrated (2.7) in the case where / = g and / is nonneg- 
ative. Extending (2.7) to the general case is then a routine application of linearity. □ 

Proof of (2.9) and the rest of Theorem 2.1. Suppose cr 2 (£) > and take / G 
1Z(Qi) with J Qi f 2 (x)dx > 0. We prove asymptotic normality for (/, /J x ), with a rate 
of convergence. To do this we adapt the proof of Corollary 2.4 of [19] (Corollary 2.1 in 
the electronically available version of [19]), to the setting of functionals of TLF point 
sets in M. d . The proof of Corollary 2.4 of [19] involves applying Stein's method to a 
graph whose vertices are sub-cubes of the unit cube with edge length proportional 
to (logA)A _1 / d and with edges between sub-cubes whenever the distance between 
sub-cubes is within twice the common cube edge length. We make the following 
trivial modifications to the proof of Corollary 2.4 of [19]. 

Let A be fixed and large. Subdivide Qf into V(A) := 3 d \p^ d sub-cubes C t A of 
volume \~ 1 pf, where p\ := a log A for some suitably large a, as in section four of 
[19]. For all 1 < i < V(X), put 

&•■= [ f(x)MPxriQ 1 ,dx)= f f( y nvnQ x ,dy). 

Then 

V(X) 

i=i 
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Note that ^ A is the analog of J^jLi l&jl °f Lemma 4.3 of [19] and furthermore, by 
the boundedness (2.1) of £, for g = 3 there exists K := K(q;f) < oo such that 

II£%<^p1 

Consider for all 1 < % < V(A) the events 

f) {^(j,A)<p A }, 

where -R^(j, A) is the radius of stabilization of £ at j G Z d . Let 

K(A) 

£a := f| Ei, 

i=i 

and note that P[^] < Ar(p A ), where r is as in Definition 2.1. 
Next, define the analog of T A in [19] by 

V(A) 

1=1 

and note that £^1^ and £^1.^ are independent whenever D 2 (C^,Cj) > 2\~ 1 l d p\. 
For 1 < % < V(X), define 

S, := (Var^'i])- 1 /^^ 

and put 

V(A) 

i=i 

As in [19] we define a dependency graph Ga := (V\,£\) for {S^}]^. The 
set Va consists of the sub-cubes Gj\ Cy(\) ano - edges {C^Cf) belong to E\ if 
D 2 (Cf, C}) < 2\-^ d p x . Next, in parallel with the proof of Corollary 2.4 of [19], we 
notice that: 

(i) V(X) := |Va| = 3 d Xp- x d , 

(ii) the maximal degree Da of Ga satisfies D\ < 5 d , 

(iii) for all 1 < i < V{\) we have \\Si\\ 3 < X (Var^'J])- 1 / 2 ^, 

(iv) Var = A), 
and 

(v) |Var[(/,pi)]-Var[(/,/i'i)]|<KA- 2 . 
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As in [19], we may use Stein's method to deduce a normal approximation result 
for S and then applying the estimates (iv) and (v) and following [19] verbatim we 
can turn this into a normal approximation result for (f,/i x ), i.e., in this way we 
obtain the desired rate (2.9) when / = 1. 

The normal approximation result for (f,/i x ), together with (2.7), implies that 
A -1 / 2 (/, fi x ) converges in distribution to a mean zero normal random variable with 
variance 0- 2 (£) Jj Q ^ d f 2 (x)dx. Given this, the convergence of the finite dimensional 
distributions in Theorem 2.1 is a standard application of the Cramer- Wold device. 
This completes the proof of Theorem 2.1. □ 

3 Stabilization of infinite input packing function- 
al 

In this section, we show that the random packing measures v\ and u' x described 
in Section 1 can each be expressed in terms of a suitably defined measure-valued 
functional £ of TLF point sets in R d x R + , of the general type considered in Section 
2, applied to a Poisson point process in space-time. Then we show that in both cases 
the appropriate choice of £ satisfies the exponential stabilization condition described 
in Definition 2.1, so that Theorem 2.1 is applicable to this choice of £. We defer 
to the next section the proof that in both cases the appropriate choice of £ satisfies 

a 2 (0 > o. 

Let us say that two points (x, t) and (y, u) in R d x R + are adjacent if (x + S) n 
(y + S) ^ 0. Given TLF X C R d x R + , let us first list the points of X in order 
of increasing time-marks using the lexicographic ordering on R d as a tie-breaker 
in the case of any pairs of points of X with equal time-marks. Then consider the 
points of X in the order of the list; let the first point in the list be accepted, and let 
each subsequent point be accepted if it is not adjacent to any previously accepted 
point of X; otherwise let it be rejected. We call this the usual rule for packing 
points of X, since it corresponds to the packing rule of Section 1 with the input 
ordering determined by time-marks. Let A(X) denote the subset of X consisting of 
all accepted points when the points of X are packed according to the usual rule. 

We consider two specific measure-valued functionals and £' on TLF point 
sets in M. d x R + , of the general type considered in Section 2, which are defined as 
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follows. For any TLF point set X C M. d x R + and bounded Borel A C R d , recall 
that A + := A x R + . Let £*(X, A) be the number of points of A(X) which lie in A + , 
and with | • | denoting Lebesgue measure, let 



e(X,A):=\S\ 



-i 



An[ |J (x + s) 

(x,t)eA(x) 



Then £* and £' are clearly translation invariant, and are bounded (i.e., satisfy (2.1)), 
since only a bounded number of solids can be packed in any fixed bounded cube. 

Recall that V\ denotes a homogeneous Poisson point process of intensity A on 
R d x R + , and V = V±. Assume V x is obtained from V by V x := \~ 1/d V. For all 
A > 0, recall the definition of £a in Section 2, and define the random measures 

• ) : = a(^A n (Q0+, •) and /if ( ■ ) := i' x {V x H (Q0+, •)■ 
Let iV| denote the total mass of /x^ , i.e. 

iVf :=Ai f(P A n[o,i]t[o, if). 

Then /i^ and /i^ are the random packing point measure and the random packing 
volume measure, respectively, corresponding to the random sequential adsorption 
process obtained by taking the spatial locations of the points of V fl Q\, in order 
of increasing time-mark, as the input sequence. Since these spatial locations are 
independent and uniformly distributed on Q\, we have the distributional equalities 

/if = u x , 4 I u' x , and iVf I iV A , (3.1) 

where the measures v\ and z/ A are given in (1.3) and the jamming number N\ is also 
given in Section I. 

We show in Lemmas 3.4 and 3.6 below that both £* and £' are exponentially 
stabilizing, and therefore we can apply Theorem 2.1 to either of these choices of £. 
To proceed with the proof of exponential stabilization, consider a partition of M. d 
into translates of the unit cube C :— Q\ — [0, \) d . It is convenient to index these 
translates as Cj, i := (ii, . . . ,id) £ with Cj := (ii, . . . , id) + C. We shall write 
C t + := Uj eZ d ||j_j|| 00 <i Cj) that is to say Cf is the union of Cj and its neighboring 
cubes. We also consider the moat ACi := C t + \ Cj. 
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We need further terminology. Given TLF X C R d x R + , and given A C R d , we 
say that fully packs the region A if every point in A + is adjacent to at least one 
point of A(X). For t > 0, we say X fully packs A by time ttf X n(R d x [0, t]) fully 
packs A. Given B C M d , we say that a finite point configuration A" C (B fl C+)+ is 
maximal or strongly saturates the cube Cj in £? if for each TLF external configuration 
y C (B \ Cj + )+, rYU}^ fully packs the region B fl Cj (the existence of maximal 
configurations is guaranteed by Lemmas 3.1 and 3.3 below). 

We shall be interested in strong saturation of Ci in B when 5 = M d or when 

= Q\. The reason for our interest is this: If we knew that there was a constant 
t < oo such that V fl (Cq" x [0, t]) strongly saturated C in M d a.s., then points 
in P with time marks exceeding r would have no bearing on the packing status of 
points in V fl (Co)+. Thus, to check stabilization of £ at it would be enough to 
replace V by the Poisson point process V fl (M d x [0, r]), and follow the stabilization 
arguments for packing with finite Poisson input (section four of [17]). While clearly 
no such constant r exists, we shall show in Lemma 3.3 that a finite random r exists. 

We say that X locally strongly saturates C{ if for each r\ C X(l (ACj) + , the point 
set (X fl (Cj)+) U ?7 fully packs Cj. The following lemma shows that local strong 
saturation implies strong saturation. 

Lemma 3.1 Suppose X C (Cj + ) + TLF and locally strongly saturates Ci. Then 
for any BCf 1 with Ci C B, X H B strongly saturates Ci in B. 

Proof. Let y C (B \ C+)+ be TLF. Let 77 := fl B + ) U y) n (AC;) + . We claim 

that 

«4((^ n b + ) u y) n (C+)+ = ^((A" n (C t )+) u 77). (3.2) 

Indeed, considering each point of (X fl (Cj)+) U 77 in the usual temporal order, we 
see that the decision on whether to accept is the same for these points whether we 
are applying the usual packing rule to (X fl B + ) U y or to (X fl (Cj)+) U 77. 

Since we assume X locally strongly saturates Ci, (X fl (Ci)+) U 77 fully packs Cj, 
and so by (3.2), (# n B+) U y fully packs Q. □ 

We will use one more auxiliary lemma. 
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Lemma 3.2 With probability 1, V has the property that for any r) C V fl (ACo) + , 
there exists T < oo such that the point set (V fl (C )+) Ur? /uZ/y pacA;s C fry time T. 

Proof. Suppose that for each rational hypercube Q contained in Co, V fl Q+ ^ 0; 
this event has probability 1. 

Take 77 C V n (AC )+. Let A. := fl (C )+) U 77). Clearly .A is finite. Let 

V be the set of x G Co such that (x, 0) does not lie adjacent to any point of A. 
Then V is open in Co (because we assume S is closed) and if it is non-empty, it 
contains a rational cube contained in Co so that V + contains a point of V n (Co)+. 
But then this point should have been accepted so there is a contradiction. Hence V 
is empty and since A is finite this shows that C is fully packed within a finite time. □ 

For i G Z d , let Tj := Ti(V) denote the time till local strong saturation, defined to 
be the smallest t G [0, 00] such that C is locally strongly saturated by the point set 
(Pn(C+) + )n(R d x [0,t]) (and set T { = 00 if no such t exists). Clearly, Tj, % G Z d , are 
identically distributed random variables depending only on Vn(C^) + . In particular, 
(Tj, i G Z d ) forms a 2-dependent random field, meaning that T { is independent of 
(Tj, \\j — iHoo > 2) for each % G Z d . We can now prove the key result that T is 
almost surely finite. 

Lemma 3.3 It is the case that P[T = 00] = 0. 

Proof. Suppose that To = 00. Then for each positive integer r there exists i] T C 

V fl (ACo) + such that (V fl (C )+) U r\ T does not fully pack C by time r. 

Assume V fl (ACo)+ is locally finite (this happens almost surely). Then V fl 
(ACo x [0, 1]) is finite so that we can take a subsequence r' — > 00 of r along which 
r] T /fl(ACo x [0, 1]) is the same for all r'. Then we can take a further subsequence r" of 
t' along which r\ T n n (AC x [0, 2]) is the same for all r". Repeating this procedure 
and using Cantor's diagonal argument, we can find a subsequence r n tending to 
infinity, and a limit set rj C (AC x R + ), such that for all fc, it is the case that 

r) Tn n (ACo x [0, A;]) = 7] n (ACo x [0, k}) (3.3) 

for all but finitely many n. 



19 



Let k > 0, and choose n to be large enough so that r n > k and such that (3.3) 
holds. Then the point set (V fl (C )+) U r] Tn does not yet fully pack C by time r n , 
and therefore (V fl (C )+) U rj does not yet fully pack C by time fc. 

Since (V fl (C )+) U r\ does not yet fully pack C by time for any k, we are 
in the complement of the event described in Lemma 3.2. Thus by that result, the 
event {T = oo} is contained in an event of probability zero, which completes the 
proof of Lemma 3.3. □ 

Using Lemma 3.3, we can now prove that £* and defined at the start of this 
section, satisfy the first part of exponential stabilization (exponential decay of the 
tail of R'). 

Lemma 3.4 There exists a positive constant K\ such that for either £ = £* or 
£ = there is a stabilization radius R' as described in Definition 2.1, satisfying 

P [R' > L] < Ki exp(-L/X 1 ), VL > 0. 

Proof. Let 5± > be a number falling below the critical probability p c for site 
percolation on Z d with neighborhood relation i = (i 1? . . . ~ j = (j 1? . . . , j^) if 
and only if ||i — J 1 1 oo ^ 1) see Grimmett [8]. 

We will apply a domination by product measures result of [9], more precisely 
Theorem 0.0 in [9]. This tells us that, for a family of {0, l}-valued random vari- 
ables indexed by lattice vertices, if we are able to show that for each given site 
the probability of seeing 1 there conditioned on the configuration outside a fixed 
size neighborhood of the site exceeds certain large enough p, then this random field 
dominates a product measure with positive density q which can be made arbitrarily 
close to 1 by appropriate choice of p. By this result, with 5i as chosen above we 
can find 5 2 > such that any 2-dependent random field (Y i: i G Z d ) with Y { taking 
values in {0, 1} and P[Yi = 1] > 1 — 5 2 for each i, this random field dominates the 
product measure with density 1 — Si. 

Using Lemma 3.3, take T* > such that P[T > T*\ < 5 2 . Then by the 
conclusion of the preceding paragraph, V can be coupled on a common probability 
space with an i.i.d. {0, l}-valued random field n i: i G Z d , so that, for all i G Z d , 

• P[7Ti = 1] > 1 - Si, 
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• we have Tj < T* whenever 7Tj = 1. 

Let us say that the cube Cj is T*-saturated if T« < T*. By Lemma 3.1, if Cj is 
T*-saturated then for any BC| d with QC B,Vf) ([C+ n £] x [0,T*]) strongly 
saturates Cj in B. 

We declare a point (x,t) G P fl (Cj)+ to be causally relevant if either 

• 7Tj = 0, 

• or 7Tj = 1 and t <T*. 

Otherwise the point x E V P\ (Cj)+ is declared causally irrelevant. 

We now argue as follows, directly adapting the oriented percolation based tech- 
nique introduced in section four of [17]. We convert the collection of points V (in 
R d x R + ) into a directed graph by providing a directed connection from (y, s) to 
(x,t) whenever |y — x| < 2d$ and s < i and, moreover, both (x,t) and (y,s) are 
causally relevant. By the causal cluster Cl[(x, t); V] of (x, t) G V we understand the 
set of all causally relevant points (y, s) of P such that there is a directed path from 
(y, s) to (x, t) (referred to as a causal chain for (x, t) in the sequel). Necessarily the 
points in the causal cluster for (x, t) have time mark at most t. 

For each (x, t) G V we define the causal cube cluster of (x, t) in ~R d by 

C1[(:M);P] :=lJ[Cy- : (C,-)+ H Cl[(x, t); P] ^ 0] 

and for each j 6 Z d we define its causal cube cluster as the union of clusters given 
by 

Cl[i;P]:= |J C1[0M);P]. (3.4) 

(x,t)e7»n(C^) + 

The significance of causal cube clusters is as follows. First, we assert that 
the packing status of a given point (x, t) is unaffected by changes to V outside 
C\[(x,t);V}. Indeed, viewing the directed connections as potential direct interac- 
tions between overlapping solids in the course of the sequential packing process, we 
can repeat the corresponding argument from Lemma 4.1 in [17], adding the extra 
observation that causally irrelevant points will not be accepted regardless of the 
outside packing configuration and hence do not have to be taken into account. Sim- 
ilarly, the packing status of the totality of points falling within distance ds of the 
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cube Ci can only be affected by the status of points falling in the causal cube cluster 
Cl[z; V]. Consequently, we see that for either £ = £* or £ = £', we can define a radius 
of stabilization by 

R' :=diam(Cl[0;P]). (3.5) 

We need to show that R' is almost surely finite with an exponentially decaying 
tail. Given L > 0, let E 1 (L) be the event that there is a 'path of zeros' from some 
site % G { — 1,0, l} d to the complement of B L , 2 _^ in the Bernoulli random field 
(7Tj,i G Z d ). More formally, Ei(L) is the event that there exists there is a sequence 
iq, i-y, i 2 , . . . , i n , such that (a) i G {—1,0, l} d , and (b) % n £ B L ^ 2 _ 2 ^, and (c) for 
j — 1, . . . , n, ij G Z d and — ij-i||oo = 1 and ir^ = 0. 

For % G Z d , let E 2 (L,i) be the event that there exists (x,t) £ V (~) (Cj) + , such 
that t < T* and there exists a causal chain for (x, t) which starts at some point of 
V \ (B L _ 2y /2) + . Define the event 

E 2 (L) := |J [E 2 (L, i):ie Z d , Q fl B L/2 ^ 0] . 

Then we assert that the event {R' > L} is contained in Ei(L) U E 2 (L). Indeed, if 
E 2 (L) does not occur, then for any causal chain for any (x,t) G V fl (Cq) + starting 
outside (B L _ 2 ^j) + , all points in the causal chain of (x,t) lying inside (B L / 2 ) + must 
have time-coordinate greater than T*; if also Ei(L) does not occur, at least one 
of these points must lie in a cube which is T*-saturated, and therefore be causally 
irrelevant, so in fact there is no causal chain for any (x,t) G V fl (Cq) + starting 
outside (B L _ 2V j) + . Hence, C\[0,V] C B L , so that R! < L. 

By the choice of 5\ and by the exponential decay of the cluster size in the 
subcritical percolation regime (see e.g. Sections 5.2 and 6.3 in Grimmett [8]), we 
have exponential decay of P[E\(L)\. That is, there is a constant K 2 such that 
P[£i(L)] < K 2 exp(-L/K 2 ) for all L. 

Since T* is fixed, we can use the methods of [17] for finite (Poisson) input pack- 
ing, in particular the argument leading to Lemma 4.2 in [17], to see that there is a 
constant K 3 such that P[E 2 (L,i)} < K 3 exp(-L/K 3 ) for all % G Z d nB L/2 . Since the 
number of such % is only 0(L d ), we see that P[E 2 (L)] also decays exponentially in 
L, and hence so does P[Ei(L)] + P[E 2 (L)]. Since the event {R! > L} is contained 
in Ei(L) U E 2 (L), the lemma is proved. □ 
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To finish checking that £* and £' satisfy the conditions for Theorem 2.1, we con- 
sider strong saturation, not only of unit cubes but of cubes of slightly less than unit 
size. Let denote the cube [— 2( 1 / d ) d , i.e. the cube of side 3^^ concen- 
tric with Qq. Let us say that Qq is locally strongly saturated by a finite point set 
X C if for every 77 C X n (Qj \ Q f )+, the point set (# n ((?<)+) U 77 fully 

packs Q^. 

Lemma 3.5 Given 5 > 0, there exist constants e > and t < 00 suc/i that for all 

(e [l-e,l], 

P[P n (Qj x [0,t ]) locally strongly saturates Q c ] > 1 - 5. (3.6) 

Proof. By Lemma 3.3, we can choose t such that P fl (Qf x [0, to]) locally strongly 
saturates Q±, with probability at least 1 — 8/2. Having chosen t in this way, we can 
then choose e, with 2d s < (1 — e) 1 ^, so that for any ( E [1 — e, 1], 

p[pn((g 1 \g c )x[o,g)^0]< 5/2. 

For C < 1 with 2d s < ( 1/d , if Pn (C# x [0, t ]) strongly saturates Q 1} and ?n ((Qi \ 
Qc) x [0, to]) is empty, then V fl (Q^ x [0,t ]) strongly saturates Qc_. Hence, the 
preceding probability estimates complete the proof. □ 

Lemma 3.6 There exists a positive constant K 4 such that for either £ = £* or 
£ = t/iere is a family of stabilization radii R(i, A) = R^(i, X), defined for A > 1 
and i G 7L d as described in Definition 2.1, which satisfy 

sup P [R(i, A) > L] < K 4 exp(-L/K 4 ). (3.7) 

Proof. First let us restrict attention to A with A 1 ^ G N. Adapting notation from 
the preceding proof, for (x,t) G Vd (Q\)+ we let C\[(x, t); V fl (Qx)+] denote the set 
of all causally relevant points (y, s) of V fl (<5a)+ such that there is a directed path 
from (y, s) to (x, £), with all points in the path lying inside (Q\)+- Then define the 
causal cube cluster in Q\ for (x, t) by 

Cl[(x, t);Pn (Qx)+] := \J[C+ n Qx : (C,-)+ fl Cl[(x, t); V n (Q A ) + ] ^ 0] 
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and and for % E 7L d by 

a[<; p n (Qa)+] = |J a[(x, *); v n (Qa)4 

Or,t)ePn(Q A nc+) + 

Define 

A) := diam(Cl[i; P n (Qa)+D, \ 1/d E N. (3.8) 

Then for i e Z d , the packing statuses of points oiVCi (C+ fl <5a) + are unaffected by 
changes to V fl (<5a)+ i n the region (Q x \ B R ^ ij x)(i))+, by the same argument as in 
the preceding proof. Here we are using the fact that A 1 ^ E Z, and that if Cj C Qa 
is T*-saturated then Cj is strongly saturated in Qa by T 3 fl (Q\ x [0,T*]) (Lemma 
3.1). Thus, R(i, A) serves as a radius of stabilization in the sense of Definition 2.1 
(for either or £')■ Moreover, Cl[i; P n (Qa)+] C Cl[i; P], and so with K x as in the 
the preceding proof we have P[R(i,\) > L] < K x exp(— L/Ki), uniformly over i,X 
with \ l ' d e N. 

Now suppose A 1 ^ ^ N. In this case, instead of dividing Q\ into cubes of side 
1, some of which would not fit exactly, we divide Q\ into cubes of side slightly less 
than 1, which do fit exactly, and repeat the above argument. 

More precisely, we modify the proof of Lemma 3.4. With 62 as in that proof, 
we use Lemma 3.5 to choose constants e > and T* < 00 (with max(2<is, 1/2) < 
(1 — s) l l d ) in such a way that for any ( E [1 — e, 1] we have 

P[P n (Q+ x [0,T*]) locally strongly saturates Q c ] > 1 - 5 2 . 

With 5 thus fixed, for all large enough A we can choose ( = ((X) G [1 — e, 1] in such 
a way that \ 1 / d /^ 1 / d is an integer. Partitioning M. d into cubes C[ of volume £, we 
can then follow the argument already given for the case \ x l d E N, using the fact that 
the each of the unit cubes i + Qi, for which we need to check conditions in Theorem 
2.1, is contained in the union of at most 2 d cubes in the partition {Cj}. □ 

4 Jamming variability, variance asymptotics 

At the end of this section, we complete the proofs of Theorems 1.1 and 1.2. First, 
we need to show that the limiting variance a 2 (S,d) is non-zero for all d and all 
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S. This is achieved by Proposition 4.1 and Lemma 4.1 below. The first of these 
results establishes that any convex S C ~R d with nonempty interior satisfies jamming 
variability (as defined in remark 6, Section 1), and the second establishes that this 
is sufficient to guarantee that cr 2 (S, d) > 0. Recall from (3.1) that we can work just 
as well with iVj" as with N\. 

Proposition 4.1 The convex body S has jamming variability. 

Proof. Given S, for all iGK d define 

\\x\\ := sup{a > : (x + aS) n aS = 0}. 

It is straightforward to verify that || • || is a norm on ~R d , using the convexity of S to 
verify the triangle inequality. For nonempty A C R d , and x G R d , write D(x, A) for 
inf { 1 1 x — y\\ : y G A}. By our earlier assumption that 2d s < 1 we have ||x|| < ||^||oo 
for all x G M. d . 

For C C M d , we shall say £ is packed if ||x — y\\ > 1 for all x G £, y G C, and 
that £ is maximally packed if it is packed and 

D(w,C)<l, VweR d . (4.1) 

We shall say £ is a periodic set if for all x G C and z <E Z d we have i + zG £. 

Let £ be a maximally packed periodic subset of R d (it is not hard to see that 
such an £ exists). Then the function x h- > £) is a continuous function on R d 
that is periodic (i.e., D(x,C) = D(x + z,£) for all x e R d ,z e Z d ). Hence the 
range of this function is the continuous image of the compact torus W d /Z d , and so 
is compact. Hence by (4.1) we have 

P := sup{D(w, £) : w G R d } < 1. 

Then for x G M d and a > 0, by scaling 

D(x, aC) = aD(a- l x, £) < af3. (4.2) 

Choose 5 > such that (3(1 + 65) < 1 - 25. For i = 1, 2, let £* := (1 + 3i<S)£. By 
(4.2) and the choice of 5 we have for all x £M. d and i — 1,2 that 

D(ar, A) < 1 - 25. (4.3) 
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Let ci denote the number of points of C in [0, l) d . Denote by Box(L) the hy- 
percube [-L/2, L/2] d . For % — 1,2, let rij(L) denote the number of points of Ci in 
Box(L — 4). Then as L — > oo, for i = 1, 2 we have 

n;(L) ~ ci(l + 35i) _d L d . (4.4) 

Let ri3(L) denote the maximum integer m such that there exists a packed subset of 
Box(L) \ Box(L — 6) with m elements. Then there is a finite constant c 2 such that 
for all L > 6 we have 

n 3 (L) < c 2 L d - 1 . (4.5) 

By (4.4) and (4.5), we can choose L such that for L > L we have 

n 3 (L) < ni (L)-n 2 (L). (4.6) 

For iGR d and r > 0, set B r (x) := {y G M. d : \\y — x\\ < r} (a ball of radius r using 
the norm || • ||). For bounded A C M. d , let T(A) denote the time of the first Poisson 
arrival in A, i.e set 

T(A) := inf {t : V D (A x {t}) ^ 0}, 

with the convention that the infimum of the empty set is oo. Fix L > L Q , and for 
i — 1,2 define the event E,- L by 

Ei := {m&x{T(B 6 (x)) :xe&n Box(L — 4)} < T (Box(L) \ U xeCi nBax(L-4)Bs(x)) } ■ 

Let % — 1 or % — 2. If y,y' are distinct points of Ci then \\y — y'\\ > 1 + 35. Hence, 
if also w G Bs(y) and w' G B$(y'), then ||w — > 1 + 5 by the triangle inequality. 
Moreover, for x G M d , by (4.3) and the triangle inequality we can find y = y(x) G Ci 
such that \\x — w\\ < 1 — 5 for all w G Bs(y). Hence, if Ei occurs then the set 
of accepted points (i.e., centroids of accepted shapes) of the infinite input packing 
process on Box(L) induced by V with arbitrary external pre-packed configuration r\ 
in R d \ Box(L), includes one point from each Bs(x),x G^fl Box(L — 4), and also 
contains no other points from Box(L — 6). 

Thus for any pre-packed configuration rj in M d \Box(L), if E\ occurs the number 
of accepted points in Box(L) is at least ni(L), and if E 2 occurs the number of 
accepted points is at most n 2 (L) + n 3 (L). Also, the probabilities P[E\\ and P[E 2 ] 
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are strictly positive and do not depend on rj. By (4.6), it follows that there is a 
constant e > 0, independent of rj, such that Var[iV"| d (Box(L))|r7] > e. Thus we have 
established the required jamming variability. □ 

Lemma 4.1 It is the case that liminfA^oo A _1 Var[iVf ] > 0. 

Proof. By Proposition 4.1, there exists L > such that inf^ VariV[[0, L] d \rf\ > 0, 
where the infimum is over all admissible rj C M. d \ [0, L] d . We consider A with 
\ 1/d /(L + 4) e N. We subdivide the cube Q\ into n(A) := A/(L + 4) d equal- 
sized sub-cubes C^a, C^a, • • • , C^a),* arising as translates of Box(L + 4) centered at 
x lt \, . . . ,x„(a),a respectively. For 1 < % < n(A), let C~ x be the translate of Box(L) 
centered at x i: \, and let M iy \ be the translate of Box(L + 2) \ Box(L) centered at 
Xj A (a 'moat' around C^ x )- 

Using terminology from Section 3, let F^a be the event that the point set V H 
(Mi ! x)+ fully packs M i: \ by time 1, and let G iy \ be the event that V C\((C^\\M^\) x 
[0, 1]) is empty. Let E it \ := F i>x n G^a- Then p := P[E it \\ satisfies p > 0, and does 
not depend on i or A. 

Observing that the events Ei^\,l < i < n(X), are independent (the cubes C{ 
are disjoint), denote the (random) set of indices for which E^a occurs by /(A) := 
{ii, ixix)}- Then E [if (A)] = pn(X). Conditional on the event E it \, the packing 
process inside C~ x has a particularly simple form - before time 1 there are no points 
in C~ x , and after that time the newly arriving solids centered in C~ x undergo the 
packing process according to the usual rules with the additional restriction that a 
solid overlapping another one packed in M it \ before time 1 is rejected. Note that 
for i G i(A), no new solids are accepted in M it \ after time 1 and, moreover, the 
acceptance times of solids accepted in M it \ before time 1 have no influence on the 
behavior of the packing process in C~ x after time 1; only their spatial locations 
matter. For a configuration rj of accepted points (only spatial locations taken into 
account) in M ijX , the process described above will be referred to as packing in C~ x 
in the presence of the pre-packed configuration rj. 

Let M\ be the sigma-algebra generated by the points of V fl (Q\ x [0, 1]), i.e. 
the Poisson arrivals up to time 1. Event E iy \ is A^A-measurable, for each i. 

By the conditional variance formula we have 



Var 



N x =E Var(iVfLM A ) + Var E^7Vf|A^ A ) 
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where we set iV| [C^Tj := H Q\,C^ X ), the number of solids packed in C^. 

Conditionally on Ai\, the packing processes after time 1 over different sub-cubes 
C~ x , i G /(A), are independent of each other and of the packing process after time 1 
in Qx \ U ieIW C~ x . Hence, 



Var 



JVf] >E E [Var[JVf [C~ x ] \ M x ]\ > E [K] inf Var7V[[0, L] d \rj\, 



where the infimum is taken over all admissible configurations rj outside [0,L] d , and 
where N[[0, L] d \rj] stands for number of solids packed in [0,L] d in the presence of 
the pre-packed configuration rj. By Proposition 4.1, this infimum is strictly positive, 
and Lemma 4. 1 follows. □ 



Proof of Theorems 1.1 and 1.2. Let £ be £* as defined in Section 3. Then Lemmas 
3.4 and 3.6 show that £ = £* satisfies the exponential stabilization conditions in 
Theorem 2.1, so it satisfies the conclusions (2.6), (2.7) and (2.8) of that result. The 
conclusion (2.8) gives us (1.1) of Theorem 1.1. Also, by putting / = g = 1 on Q\ 
and using (2.7), we obtain the variance convergence A _1 VariVA — > o 2 asserted in 
Theorem 1.1. By Lemma 4.1, we may therefore deduce that a 2 > 0. Hence we may 
apply the last part of Theorem 2.1 to obtain the rest of the conclusions in Theorem 
1.2 as they pertain to u x ; also the conclusion (2.9) of Theorem 2.1 gives us (1.2). 

To get the same results for i/, we argue similarly with £ = We need to 
check that the limiting means and variances are the same, i.e. M£') = /-*(£*) and 
<7 2 (£') = (T 2 (£*). To see this, note that if / = 1 on Qf, then (/,/xf) = so 
application of (2.6) to this choice of / yields 

MO = lim A- X E [</,/4'>] = lim A _1 E [(/,//>] = MO 
and a similar argument using (2.7) shows that cr 2 (£') = cx 2 (0- □ 
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